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^ ■ Abstract: We compute the prediction of the Nambu-Goto effective string model for a 

^^ ■ rectangular Wilson loop up to three loops. This is done through the use of an operatorial, 

^r^ ■ first order formulation and of the open string analogues of boundary states. This result is 

''^ . interesting since there are universality theorems stating that the predictions up to three 

t^ \ loops are common to all effective string models. To test the effective string prediction, we 

use the Montecarlo evaluation, in the 3d Ising gauge model, of an observable (the ratio 
of two Wilson loops with the same perimeter) for which boundary effects are relatively 
small. Our simulation attains a level of precision which is sufficient to test the two-loop 
^ ■ correction. The three-loop correction seems to go in the right direction, but is actually yet 

beyond the reach of our simulation, since its effect is comparable with the statistical errors 
of the latter. 
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1. Introduction 

The proposal that the strong-coupling dynamics of gauge theories could be captured by 
an effective string theory is a long-standing one [1]-[10]; the fluctuating long string, in 
this approach, is the color flux tube joining coloured charges. Perhaps the most direct 
prediction that can be extracted from the effective string is the shape of the potential 
V{R) between two external sources (two static quarks) at distance R. From the gauge 
theory point of view, if we have a rectangular Wilson loop W{L, R) of sides L and R, the 
inter-quark potential is 

V{R) = - lim ^logW{L,R) . (1.1) 

In the confining phase, the area law for the Wilson loop corresponds to a linear potential 

V{R) = aR+... . (1.2) 

In a string interpretation, the area term aLR in the exponent of the Wilson loop is the 
classical action of an effective string model. Another way of measuring the static qq poten- 
tial is to consider the correlator of two Polyakov loops at distance i? in a gauge theory at 



finite temperature oc 1/L; in the strong-coupling regime, the string world-sheet describing 
this situation is a cylinder. 

The simplest and most obvious choice for the effective string is the Nambu-Goto one 
[11, 12], where the action is the induced area of the string world-sheet, see eq. (2.1) below, 
and the degrees of freedom are the D = d — 2 transverse fields. The action is non linear and 
contains, besides the kinetic term, a series of higher derivative interactions. Quantization 
is implemented by functional integration and can be carried out perturbatively. This 
approach was used in [13] to compute, up to two loop order, the string vacuum amplitude 
on world-sheets with disk, cylinder or torus topology (the first two being relevant for Wilson 
loops and Polyakov loop correlators respectively). 

On general grounds, one expects any effective theory for a long string to contain as 
degrees of freedom the transverse fields, representing the Goldstone modes for the transla- 
tional invariance broken by the string configuration. The "one loop" quantum corrections 
are then given by the functional determinants of the kinetic operators for the transverse 
modes and are the same for all effective models. They lead to the correction to the static 
potential known as "Liischer term" [9, 10]: 

V{R) = aR-^ + ... . (1.3) 

Over the years, computer simulations of Polyakov loop correlators and Wilson loops have 
shown with increasing evidence the presence of this correction (for an up to date review 
see [14]), confirming the soundness of the effective string approach. 

The various possible effective theories are distinguished by the form of their interaction 
terms, which in turn affect the perturbative corrections to the amplitudes, starting at two 
loops. The interaction terms cannot, however, be completely arbitrary. In [15] it was shown 
that the requirement that the effective string cylinder partition function (corresponding to 
a Polyakov loop correlator) can be re-expressed in terms of propagating closed string states 
fixes the coefficients of the first higher derivative terms. Moreover, if the effective string 
is supposed to describe the low-energy regime of a consistent quantum field theory, the 
Poincare symmetry of the latter must be still realized. Since the configuration around which 
the long string fluctuates spontaneously breaks, in particular, the Lorentz transformations 
mixing longitudinal and transverse coordinates, these must be realized non-linear ly in the 
effective action. This requirement is present in the analysis of [15], see [16], but was made 
more explicit and pursued further in [17]- [19], where it is shown that the coefficients of 
the quartic and sextic derivative interaction terms are almost completely fixed by this 
constraint, both for closed and open strings, and coincide with those of the NG model. 
This "universality" of the lowest interaction terms suggests that Montecarlo simulations 
should be compatible with NG prediction also beyond the one-loop level. 

Testing the effective string predictions beyond the one-loop terms, for instance by look- 
ing at subleading terms in the 1/R expansion of the effective potential, is a very serious 
challenge. It is a challenge from the computational point of view, because it requires to 
carry out simulations so precise as to distinguish deviations of the results from the one-loop 
predictions. However, thanks to the development of new algorithms and the surge in avail- 
able computer power, such a precision is basically within reach. In particular, in the case 



of the three dimensional gauge Ising model in which, exploiting duality, high precision can 
be reached at a reasonable computational cost, universality was unambiguously confirmed 
at two loop level both in the case of the torus geometry [20, 21, 22] and in the case of the 
cylinder geometry [23, 24]. Testing universality at three loop level turned out to be more 
difficult. The only existing tests [25] were performed in very asymmetric geometries were 
string effects are magnified but systematic errors, due for instance to the vicinity of the 
deconfinement transition, are more difficult to control. For a survey of most recent results 
for other gauge groups or other observables see [14]). Until now no attempt to go beyond 
the one loop results was performed in the case of the disk geometry. Filling this gap is one 
of the main goals of the present paper. 

Also on the theoretical side, extracting the higher loop effects of the effective string 
model is not an easy task. The diagrammatic evaluation of the effective string amplitudes 
carried out in [13] was pushed up to three loops, and extended to a generic effective string 
action to this order in [17], for the cylinder partition function only. 

In this work, we derive higher loop corrections in the Nambu-Goto disk amplitude 
corresponding to the Wilson loop. Our method, which allows in fact to obtain the exact 
amplitude that resums the loop expansion, is based on the the first-order reformulation of 
the Nambu-Goto string model; it has already been successfully applied to the cylinder [26] 
and the torus [27, 28] partition functions. We will see that the Wilson loop case is more 
delicate, but we propose an exact expression which reproduces, at two-loop, the result of 
[13]: this represents a very non-trivial check and we trust therefore our results also to 
higher loops. 

The universality of the first terms in the derivative expansion of the effective action 
explains why the NG model represents a good approximation of the correct effective string 
theory for the flux tube; in turn, testing NG to the order to which its effects are universal 
would represent a stringent test of the whole effective string approach. 

Still, the NG model can not be the right one, and deviations from its predictions, when 
precisely identified, may help in the quest for the correct effective string. In particular, the 
Lorentz invariance should be preserved also at the quantum level. The quantum consistency 
of the NG theory is usually investigated by using a first-order re- formulation, which involves 
an independent metric on the world-sheet but is quadratic in the string fields: employing the 
Weyl invariance to reach the so-called " conformal gauge" , the world-sheet metric decouples 
and one is left with a quadratic action, to be supplemented with the Virasoro constraints. 
This is the standard approach described in string theory textbooks, see for instance [29, 30]. 
The residual conformal invariance generated by the modes of the Virasoro constraints can 
be exploited to fix a light-cone (physical) gauge where only the D transverse fields are 
independent; it can also be enforced, in the "covariant quantization" , by taking into account 
the effects of the ghost/anti-ghost system that arise from the Jacobian to fix the conformal 
gauge. In the light-cone gauge, the closure of the target-space Lorentz algebra at the 
quantum level requires D = 24; however, when the string fluctuates around a long string 
configuration with length scale R, this anomaly is suppressed in the large R limit [31]. In 
covariant quantization, the central charge Cgh = —26 of the ghost/anti-ghost system has to 
be compensated by that of D -|- 2 = 26 bosonic fields. 



Polyakov's approach [32] of functionally integrating over the intrinsic metric shows that 
for D ^ 2Q the scale of this metric does not decouple and represents an additional degree of 
freedom that has to be taken into account to ensure quantum consistency. This world-sheet 
boson is usually referred to as the Liouville field since its action is of the Liouville type. If 
the effective string theory must contain just the D transverse d.o.f., and realize also at the 
quantum level the Lorentz algebra, none of the standard quantizations can thus be correct. 
An interesting proposal was put forward by Polchinski and Strominger in [33] ; to reconcile 
the requirement of critical central charge with having just the D transverse degrees of 
freedom, they proposed a Polyakov model where the Liouville field is re-expressed in terms 
of the induced metric, i.e., in terms of the transverse fields. The action becomes non-local, 
but admits nevertheless a sensible expansion around long string configurations. It has been 
shown that the first terms in the derivative expansion of this model agree with the NG 
ones, in accordance with the universality arguments. It seems however difficult to compute 
higher loop corrections to the amplitudes in this model. 

Recently, some particular, supersymmetric, gauge theories have been finally given an 
explicit realization in terms of strings propagating on a curved manifold with an extra 
dimension, in the so-called AdS/CFT duality [34, 35, 36]; the extra dimension can be 
interpreted in terms of the energy scale [37] and as such it does not represent a spurious 
degree of freedom. From the point of view of the effective action for the transverse fields, 
the interaction terms would take into account the curved geometry. This is of course a 
very intriguing possibility. 

There is an extra issue in devising the correct effective string action, namely the pos- 
sible presence of boundary terms. Lattice simulations have made it clear that the Wilson 
loop expectation value displays, beside the leading area term, a perimeter term with an 
independent coupling: 

W{L, R) ~ exp {-aLR - ^{L + R) + . . .) . (1.4) 

In this paper, we will study a particular observable, corresponding to the ratio of two 
different rectangular Wilson loops, such that the leading perimeter dependence described 
in eq. (1.4) cancels [38]. However, it is natural to expect that the perimeter term is 
only the classical value of some boundary component of the effective action which also 
yields other, subleading, corrections to the amplitudes. In [15, 18, 19] the first few possible 
boundary terms in a derivative expansion of the effective action have been investigated 
for the Polyakov loop geometry, finding that they can lead to corrections of the order of 
1/i?^ {R being the distance between the Polyakov loops); to our knowledge, no explicit 
evaluation of the effect of boundary terms in the Wilson loop geometry is available. From 
the theoretical point of view, in [39]- [43] by considering Polyakov quantization in presence 
of boundaries, it was found that boundary terms in the world-sheet actions are needed, 
and in particular a term proportional to the induced length of the boundary, whose leading 
contribution to the Wilson loop amplitude would indeed be a perimeter term. It would be 
nice to investigate the modifications of the Nambu-Goto calculations in presence of such a 
term, but we leave this to future work. 



In comparing the outcomes of our simulations to the theoretical predictions of the 
Nambu-Goto effective model, we will have to take into account the likely presence of bound- 
ary terms, of which no explicit evaluation is available; this increases the uncertainties in 
the comparison. 

As mentioned above the main goal of this paper was to compare theoretical predictions 
and Montecarlo simulations beyond one loop in the case of the Wilson loop geometry. We 
shall show that, using duality based algorithms, it is possible to study Wilson loops large 
enough to make boundary effects negligible and at the same time to reach a precision high 
enough to disentangle effective string corrections beyond the one loop level. As in the cases 
of the torus and cylinder geometry, we shall be able to confirm universality at the two loops 
level. The third loop corrections seem to go in the right direction, but their magnitude is 
of the same order as the statistical errors, so that they remain beyond the reach of present 
simulations. 

2. Two- loop effective string prediction 

In this section we summarize the two-loop prediction for the Wilson loop amplitude ob- 
tained in the Nambu-Goto effective string model. As discussed in the introduction, the 
predictions up to three loops are actually universal for all effective string models. 

In the Nambu-Goto approach [11, 12] the action for the fluctuating surface spanned 
by the flux tube is simply its induced area in the d-dimensional target space: 

(f^y^detg , ga(5 = -^t^-^t^Gmn ■ (2.1) 

Here a is the string tension, the surface is parametrized by proper coordinates ^", and 
X^{^) (M = 1, . . . , D + 2) describe the target space position of a point specified by ^. 
For us the target space metric Gmn will always be the flat one. Our analysis aims at 
comparing effective string predictions with lattice simulations of gauge theories, which 
are always performed with an Euclidean time. As a consequence, we use the Euclidean 
signature on the target space as well as on the world-sheet; this is what is commonly done 
in the literature about the effective string approach 

Consider a rectangular Wilson loop of sides L and R in the x^,x^ plane. The direction 
x^ is interpreted as the Euclidean time, so that in the limit L — t- oo one can extract the static 
quark- antiquark potential via eq. (1.1) and compare it directly with lattice simulations. 

In the so-called "static" gauge the proper coordinates of the string are identified with 
x^ and x"^ already at the classical level. The partition function for the Wilson loop surface 
is then obtained by functional integration over the transverse fields X(x^,x^): 

W{L, R)= I DX exp | -^ /" dx^ f dx^ (l + {diXf + {d2Xf + {diX A d2Xf^ ' 1 



/ DX exp I -a / dx^ 



R 

dx^ 



1 + -{diXf + -{d2Xf + interactions 



(2.2) 



with the transverse fields X vanishing on the Wilson loop perimeter. To this NG action, 
terms living on the boundary of the domain can (and must) be added, and their effect is 
expected to be important, especially for small loops. In this paper, however, we will not 
deal explicitly with the quantum corrections induced by such terms. Eq. 2.2 describes a 
statistical sum over surfaces bounded by the Wilson loop. 
By construction, we must have 

W{L,R) = W{R,L) . (2.3) 

Expanding the square root as in the second line of eq. (2.2), the constant term leads 
to the well-known classical area law 

W^i{R, L) = e-^RL+P(R+L)+k _ ^2.4) 

Here on top of the area law we have included a perimeter term, whose strength is parametrized 
by p, which arises from boundary terms in the action. The overall normalization, which is 
out of control, is parametrized by k. 

From the functional integration over the fields X, with the bulk action given in the 
second line of eq. (2.2), one gets the determinant of their kinetic operator multiplied by 
perturbative corrections, starting at two loops. The result has the structure 

,,„,,, . ,-..««..,«. (*!)-*{,, ^ , 1^ ,, ,j , ,,5) 



where, for ease of notation, we have introduced 

A = RL , u=- . (2.6) 

R 

In eq. (2.5), the term {ri(\u)/\/R) ''^, where by rj we denote Dedekind's function (see 
Appendix A for relevant notation and definitions) is the functional determinant; it is in- 
variant under i? -H- L, i.e., under u -H- 1/u, as it follows from the modular transformation 
properties of the rj function. The corresponding contribution to the interquark potential of 
eq. (1.1) is the celebrated Liischer term Dti/{2AR) appearing in eq. (1.3) that represents 
the leading correction to the linear confining potential [9, 10]. 

With the impressive enhancement in precision of Montecarlo simulations in the last 
years and, above all, in view of the universality theorems discussed in the introduction it 
becomes important to evaluate and test also the subleading terms in eq. (2.5). The terms 
Ca{}u) represent the corrections at order a in the loop expansion parameter a A; they must 
be invariant under u -H- 1/n. 

The computation of the two-loop^ correction £2 was carried out almost 30 years ago 
by Dietz and Filk [13]. However, the two-loop result of [13] appears to be incompatible 
with the data obtained with various simulation techniques, see [44] for an account of recent 



^In fact, [13] considers the three cases where the string world-sheet is a a disk, a cyhnder or a torus, 
relevant for the effective description of, respectively, Wilson loops, Polyakov loop correlators and interfaces 
in a compact target space. 



attempts in this direction. This fact prompted us to a critical re-examination of the Dietz 
and Filk calculation. We found a mistake^ in their eq. (3.3). The correct result is 



C.(u) = (i 



2D u^Ei {\u) - ^^ — ^^2 {iu) E2 ii/u) 



(2.7) 



where £"2 and £'4 denote the second and fourth order Eisenstein series, see Appendix A. 
We shall show in the next section that eq. (2.7) turns out to be compatible with the 
Montecarlo results^, in accordance with its universal character. 

3. New numerical results 

In this section we compare the two-loops effective string prediction for Wilson loops with 
a set of numerical data extracted from a Montecarlo simulation of a suitable observable in 
the three-dimensional Z2 gauge model. 

3.1 Definition of the observable in the 3d Z2 gauge model 

The 3d Z2 gauge model on a cubic lattice is defined by the partition function 



ZaiP) = 2^ exp{-/35,} . (3.1) 

The action Sg is a sum over all the plaquettes of the lattice, 

'S'c, = - 2_^ (Tn , o-n = (yi^ai^ai^ai^ , (3.2) 

n 

where ai are Ising variables located on the links of the lattice. This gauge model is equiv- 
alent to the standard Ising model by the Kramers- Wannier duality transformation 

Zg{(3) cc ZsW) , (3.3) 

^ = _ilog[tanh(/3)] , (3.4) 

where Zg is the partition function of the Ising model in the dual lattice 



E 

=±1 



Z,(/3) = > ^ eM-f3Hiis)} . (3.5) 

Here 



Hi{s) = -y~] J/ij)SiSj (3.6) 






^The term proportional to A2 in that equation should read D{{3) + (4) +2(1) +4(2)) instead of -D((3) + 
(4) + 2(1) — 4(2)). This mistake was not easy to spot because its contribution vanishes in the cylinder 
and torus geometries; in these cases, in fact, the two loop prediction of Dietz and Filk is compatible with 
the Montecarlo simulations [20, 21, 23]. Moreover the term affected by this error is symmetric under the 
exchange R -^ L and thus can not be excluded using this symmetry requirement. 

^Note however that, as said before, the presence of boundary terms in the string action leads, besides 
the classical perimeter term, to further quantum corrections to eq. (2.5), which we are not taking into 
account in our theoretical predictions. 
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Figure 1: The observable R^{L,R,n) is given by the ratio of the Wilson loops W{L,R) 
and W{L + n,R — n). 



with i and j denoting the nodes of the dual lattice and the sum is over all the links {ij) 
of the dual lattice. Here the couplings Jij are fixed to the value +1 for all the links. This 
define a one-to-one mapping between the free energy densities in the thermodynamic limit. 
The dual of the Wilson loop expectation value , in which we are interested, is given by 



W{R,L) 



Zs,w0) 

ZsCP) 



(3.7) 



where in Zg^y/ all the couplings of the links that intersect the surface S enclosed by the 
Wilson loop take the value Ji^ij) = — 1. 

The form of the functional dependence of the Wilson loop on the sides i?, L is given 
in eq. (2.5), where we have to fix D = d — 2 = 1. The values of the parameters a,p and k 
appearing in the classical part of this expression are not predicted by the effective string 
theory. However p and k can be eliminated by considering ratios of Wilson loops with the 
same perimeter [38], as the ones of fig. 1. We will therefore consider the observable 



Rv/{L,R,n) 



W{L,R) 



W{L + n,R-n) 



(3.8) 



3.2 Algorithm and simulation settings 

In general, computing large Wilson loops is a difficult task because the Wilson loop expec- 
tation value decreases exponentially with the area, so the relative error increases exponen- 
tially. This is true also for the gauge Ising model but in this case, if one evaluates the dual 
observable eq.(3.7), the problem can be overcome. Our approach is essentially the same 
described in great detail in [23], adapted to the observable 3.8. In order to compute eq. 



(3.8) we factorize the ratio as follows: 

W{L,R) W{L + l,R) W{L + n,R) W{L + n,R-n + l) 



R^{L,R,n) 



W{L + 1, R) W{L + 2,R)"' W{L + n,R-l)"' W{L + n,R-n) 

(3.9) 

Every ratio of eq. (3.9) is then simulated separately, using a hierarchical algorithm as 
explained in [23]. 

In order to compare the effective string theory prediction with the numerical data we 
need a very precise estimate of the zero temperature string tension. We choose to simulate 
the model at coupling /3 = 0.75180, for which the string tension is known with very high 
precision [45] to be cr = 0.0105241(15). We performed our simulations on a lattice of size 
64 X 120 X 120 and measured only Wilson loops lying along the directions of size 120. The 
lattice sizes were chosen large enough to avoid finite size effects; in particular, we carefully 
tested that a lattice size of 64 lattice spacings in the transverse direction was enough to 
avoid finite size corrections even for the largest loops that we studied. 

3.3 Comparing the results to the two-loop predictions 

We extracted from the simulation the ratio of expectation values 



n' W{L + n,^-n) 



R^iL, -, n) = „^,^ ; 'I' ^ (3.10) 



for different values of the asymmetry parameter u = L/R > 1 and of n. The numerical 
results of the simulation are reported in Table (1). The leading classical behaviour of the 
observable R„{L,L/u,n) is given by 

R„{L,-,n) = exp|-CTnfn + L(l- -)j I. (3.11) 

In order to isolate the quantum corrections in which we are interested, exploiting the fact 
that a is known with very high precision, we define the new observable 

R'„{L, -, n) = R^{L, -, n) - exp | -an (n + L{l--))\ . (3.12) 

u u I \ u J } 

We plot in fig. (2) our data for R1„ against the quantum string corrections up to two 
loop. Looking at these figures we see that a few interesting facts emerge. 

1. The deviation from the classical behavior R^ = is immediately visible from the 
plots: string corrections do definitely affect the Wilson loop ratio. 

2. The deviation of the data is systematically larger than the one loop prediction (which 
is indicated by the full line in fig. 2). For large enough values of R, the smallest of the 
two sides of the Wilson loop, the data converge toward the two loop string correction 
(the dashed line in fig. 2) and for the largest values of L we considered {R > 20) they 
agree with it within the errors. 
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Figure 2: Plot of R^{L,L/u,n) for different values of u and n against the quantum string 
corrections, up to two loop. Note that the classical string prediction is E1„ = 0, see eq.s 
(3.11,3.12). 



3. This agreement becomes worse as R decreases. This also happens for values of R 
much larger than the flux tube thickness and thus cannot be due to the breaking of 
the effective string picture. The most probable reason is the presence of subleading 
effects due to boundary terms in the effective string action. 

The first point is a well established result: effective string corrections for the Wilson loop 
were observed for the first time in the gauge Ising model in [38] and later also in several 
other LGTs and in particular in the - physically relevant - four dimensional SU(3) LGT 
[46]. The second point, instead, is a new result in the Wilson loop geometry. 

Let us address in more detail the third point. At present there is no explicit evaluation 
for the quantum effects due to boundary terms in the Wilson loop geometry. This is a 
very interesting and non trivial issue which we plan to address in future work. However, 
extrapolating to the Wilson case the existing results [19, 47] for the cylinder geometry we 
expect that these corrections should decrease at least as 1/R'^. This behaviour seems indeed 
to fit rather well the deviations observed in fig. 2. Extrapolating the fits to larger values of 
R one can show that for R > 20 boundary corrections are fully negligible. This statement 



10 



holds true assuming for the boundary corrections any possible decreasing behaviour of the 
type 1/i?" with n > 4, assuming any value, compatible with the data at low R, for the 
amplitude of the correction and, most importantly, making no assumptions on the effective 
string corrections, except for those implied by the universality theorems. 

This makes us confident that our results for the largest values of R can be used for an 
unbiased test of our theoretical predictions. In this respect it is interesting to observe that 
in this region our data suggest that the universal three loop correction should be small or 
should have the opposite sign with respect to the two loop one. In order to clarify this 
point, however, we need to compute this correction. This is the main goal of the following 
section. 
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Table 1: Results for R„{L, R, n) as a function of the long side of the loop, for various fixed 
values of the ratio u = L/R and of n. 



4. Beyond two-loops: operatorial approach 

Having argued that it is desirable to compute the three loop correction to the Wilson loop 
amplitude in the NG model, we proceed now to do so by resorting to an approach based 
on the first order re-formulation of the model itself. This approach has been used in [26] 
and [27] for the cylinder and torus partition functions respectively. The reasons to take 
this route are two-fold. On the one hand, carrying out the derivation of the three loop 
effect in the physical gauge by extending the two-loop computations of [13] is a daunting 
task. On the other hand, the application of the operatorial, first-order formalism to the 
Wilson loop disk partition function is not difficult but also not trivial and, as we will see, 
relies on the use of the rather non-standard (but perfectly sensible) concept of open string 
boundary states introduced in [48]. The way in which the loop expansion of the Wilson 
loop amplitude emerges in this formalism is quite interesting. 

We use the first order formulation, requiring the introduction of an independent world- 
sheet metric ^ab- This metric can then be gauge- fixed exploiting reparametrization invari- 
ance, and the open string action in the conformal gauge ^ reads simply 



S = <j de df \{d^iX^'f + (%X 







^M\ 



+ s. 



h. 



(4.1) 



^In the conformal gauge the world-sheet metric is of the form 7^/3 = e'^Sa/s, and corresponds to a CFT 
of central charge Cgh. — —26. The scale factor e"^ decouples at the classical level, but this property persists 
at the quantum level only if the anomaly parametrized by the total central charge c — D — 24: vanishes. We 
will nevertheless proceed in the case of general D, according to the discussion in the introduction. 
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World-sheet 



e , 



Target space 

open string with b.c.'s 
as in cq.s {4.2-4.4) 



TT 




Rcadsorbcd 



Figure 3: A surface bordered by the rectangular Wilson loop, such as the one drawn on 
the left, can be spanned by an open string with end-points attached to the sides at x^ and 
X fixed but free to move in the x^ direction. This happens if such a string is emitted from 
the vacuum at x^ = and re-adsorbed at x^ = L. 



where ^■^ G [0, vr] parametrizes the spatial extension of the string and ^^ its proper (eu- 
clidean) time evolution. The fields X^^(^^,^^), with M = 1, . . . ,D + 2, describe the em- 
bedding of the string world-sheet in the euclidean target space and form the 2-dimensional 
CFT of D -|- 2 free bosons. The term Sg^. in sq. (4.1) is the action for the ghost and 
anti-ghost fields that arise from the Jacobian for fixing the conformal gauge. We do not 
really need here its explicit expression, see [29] or [30] for reviews. 

We can treat this theory in an "Hamiltonian" way, by selecting the coordinate ^^ as a 
(radial, euclidean) "time" coordinate. At fixed ^^, the fields X describe the embedding 
of the string in the target space; evolving in ^^ the string sweeps out a surface. 



4.1 The Dirichlet string 

We argue that such an operatorial description^ is possible also with the boundary conditions 
corresponding to a rectangular Wilson loop, as described in Fig. 3. We arbitrarily select 
one of the directions along the Wilson loop, say x^, and regard the sides of the loop along 
which x^ varies, which sit at distance R in the x^ direction, as the quark and anti-quark 
lines. Consider a open string whose end-points are attached to these two lines and are free 
to move in the x^ direction. If such a string is emitted from the vacuum at the Wilson 
loop side placed at x^ = and re-adsorbed at the side placed at x^ = L then it spans a 
surface bordered by the Wilson loop. 

In this description the directions x^ and x^ play a very asymmetric role; we must 
however ensure that the Wilson loop partition function we compute be invariant under the 
exchange of L and R. To make this picture concrete, one has to quantize the open string 



^Our operatorial approach differs from the outset from the path-integral approach to the Polyakov model 
applied to the Wilson loop topology in [43, 51]. In that case, Dirichlet boundary conditions (modified so as 
to take into account reparametrization invariance along the boundary) are imposed along the entire Wilson 
loop contour. 
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with the boundary conditions just discussed, and construct then in its Hilbert space the 
states which describe its emission from the vacuum. Such states represent the open string 
analogue of the boundary states that describe the insertion of a boundary in the closed 
string world-sheet, i.e., the emission of closed strings from the vacuum (or from D-branes) 
[49], see for instance [50] for a review. 

Let us start by imposing on the string the boundary conditions discussed above. These 
are of Neumann typs conditions at both ends in the x^ direction: 

d^2X\e,f)\^2^,^^ = 0, (4.2) 

and of Dirichlet type at both ends in the direction x^: 

X\e,0) = 0, X\e,7r) = R, (4.3) 

and in the transverse ones: 

x{e,o) = x{e,7T) = 6. (4.4) 

With these boundary conditions, the string fields admit the following expansions: 

^He.e) = !^e - ,/|f: (4^-"^' + <^e<') si„<^ , (4.5) 



ra=l 

00 



Xie,e) = -.[^y f^e-<^ + ^ie^^l sinne^ . 



2 ^ \y/n ^n 

Canonical quantization leads to the following commutation relations among the oscil- 



lators a„ 



M. 



(4.6) 



The modes of the stress-energy tensor, usually denoted as Lm,, generate the residual 
conformal transformations of the world-sheet. In particular, Lq generates the world-sheet 
dilations and corresponds to the Hamiltonian derived from the action eq. (4.1). It receives 

+■'"0 



contributions from the bosons and the ghost system: Lq = -^n + -^0 ' ^'^'^ ^^ have 



n=l 

where Nn = '^Zm '^-n ' ^n i^ the occupation number for the oscillators a^ , and d/2A is the 
((" function regularized^) normal ordering constant. The term -R^ represents the energy 



^In [52] some words of caution were raised regarding the use of ("-function rcgularization for the string 
with the present boundary conditions. However, in [26] it was shown that in the Polyakov loop correlator 
geometry, where open strings have very similar boundary conditions, the results of this procedure are 
perfectly consistent. 
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due to the stretching of the string between the two opposite sides of the Wilson loop. For 
the ghost system we have, see for instance [29], 

Lq = non-zero modes H . (4-8) 

In the computation of partition functions, the net effect of the non-zero mode ghost modes 
(which have the opposite statistics with respect to ordinary bosonic fields) is to cancel 
the non-zero-mode contribution of two bosonic fields, reconciling this treatment with the 
static gauge description where only the transverse fields are present to begin with. This 
happens when the operatorial approach is utilized to derive the effective string predictions 
for Polyakov loop correlators [26] and for interfaces [27] and this will be the case also for 
the Wilson loop. 

4.2 The Wilson loop amplitude 

We have implemented the boundary conditions along the two sides in the direction x^ of the 
Wilson loop in the definition of the Dirichlet string. We have now to enforce the boundary 
conditions at the two sides in the direction x^ as operatorial condition on suitable states in 
the Hilbert space of this string. We proceed by constructing "open string boundary states" 
[48], which we denote as |-Bop,0) and \Bop,L), such that at proper time ^-"^ = we have 

X\0,f)\Bop,0)=0, X\0,f)\Bop,L)=L\Bop,L) , (4.9) 

X\0,e)\Bop,0) = ^e\Bop,0) , X\0,e)\Bop,L) = ^e\Bop,L) , 

X{0,e)\Bop, 0) = , X{0,e)\B,p, L)=6. 

The conditions on |i?op) 0) mean that, when applied to this state, the string fields X^ {£^'^) 
describe an open string coinciding with the side of the Wilson loop at x^ = 0. Analogous 
is the meaning of the conditions on | Bop , L) . 
We propose that the following expression: 

/•oo ^^ 
W{L,R)= -(i?op,0|e-2-*^«|i?op,L) , (4.10) 

Jo ^ 

which we will call "Wilson loop amplitude", be proportional to the Wilson loop free energy 
W{L, R). This expression is obtained by integrating over the real parameter t the quantum 
mechanical kernel describing the propagation of the string from the state |i?op)0) to the 
state |i?op,0) in a proper Euclidean time time 27rt. If oj were equal to zero, eq. (4.10) 
would contain the open string propagator Lq (written via the Schwinger parametrization), 
sandwiched between two boundary states. This is the same kind of expression used to re- 
formulate the open string cylinder partition function as the tree level propagation of a closed 
string between two boundary states [49, 50] that was applied to the effective description of 
Polyakov loop correlators in [26]. 

However, in order to obtain an expression invariant under the exchange L -H- i? we 
will be forced to choose a non-zero value of co, see later eq. (4.22). We do not have 
an a priori understanding of this choice; if we accepted, however, this choice not only 
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guarantees the L -H^ R symmetry, but reproduces the results of the static gauge approach 
at the classical, one-loop and two-loop level, which represents a very non trivial check. This 
a posteriori confirmation makes us confident that we can use eq. (4.10) to derive higher 
loop contributions; this can be done with a very limited effort, compared to the static 
gauge approach. Indeed, W(L,i?) can be expressed as an infinite sum of Bessel function 
contributions; moreover, it can be rather straightforwardly expanded in inverse powers of 
l/{aLR), i.e., in sigma-model loops of the NG formulation. To see this, we now consider 
the various ingredients of the amplitude eq. (4.10). 

Open string boundary states The open string boundary states have the form 

\Bop,0) = J\f\Bop,0)(^0)\Bop,0)(^n.z.)\Bop,0){gh.) ; (4-11) 

of, course, the same decomposition applies to j-Bop)-^)- Here, AA is a normalization which 
will be irrelevant for our purposes, while the subscripts refer to the zero-mode, non-zero- 
mode and ghost sectors of the open string Hilbert space. Let us analyze in turn these 
various components. From the expansion eq. (4.5) of the string fields it follows that the 
conditions eq. (4.9) in the zero-mode sector only imply that the c.o.m. position operator 
x^ assumes the values and L on \Bop,0) and \Bop,L) respectively. The zero-mode part 
of these states is therefore given by 

\Bop, 0)(o) = S{x')\0) = ^l dqi\qi) , (4.12) 

\B,p,L)i^^^=5{x^-Lm = ^jdq[e-'<%[) , (4.13) 

where |0) is the vacuum in the zero-mode sector, on which all components p of the 
momentum operator vanish, while \qi) is an eigenstate of pi with eigenvalue qi (the other 
components of the momentum vanish). 

The non-zero mode sector of the Hilbert space corresponds to a collection of harmonic 
oscillators. The conditions eq. (4.9) imply, through the mode expansions eq. (4.5), that, 
for any n, 

{a\ + a]}j l^op, 0)(n.^.) = , [al- a^t j |5op, ^)(n.z.) = , (^a„ - 4) l^op, ^){n.z.) = ; 

(4.14) 
exactly the same holds for |i?op) -^)n.z.- These conditions are easily solved^ and we have 

|i?op,0)n.z. = \Bop,L)^,,, = (g)eK-(""')'+(""')'+(^")')|0)(^) , (4.16) 

n=l 



^For an harmonic oscillator of mass m and frequency oj, such type of conditions correspond to the 
vanishing of either the momentum p = y -^^^(a + a^) or the position x — ^\ 2^iS'- ^ '^^) operator. 
Momentum or position eigenstates can be expressed in terms of the vacuum as 



\X = X) = — ] e^^ ' V R 2 R 0) , \P = P) = \ C 2^ ' ' V mh^f Ihojm Q) 



(4.15) 
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where by |0)(„-, we denote the tensor vacua for the n-th oscillation modes in all directions, 
and the overall normalization was already included in eq. (4.11). 

To avoid excessive technicality, we do not write explicitly the form of the boundary 
state in the ghost sector 

The amplitude We can now compute the matrix element appearing in the definition 
eq. (4.10) of the Wilson amplitude. From the explicit expression of the Virasoro generator 
Lq, eq. (4.7), in the zero-mode sector we find 

(o)\-yop,U|e '-'\Bop,L}(^Q) - / -T—^r—{qi\e V /e ''i jg'i) 

./ J.TT J.TT (4.17) 

27rcr -f/2-i a L^ 

e 4t , 

t 

where we used the orthogonality of momentum eigenstates: {qi\q'i) = 2Ti5{qi — q'^) and 
performed the remaining Gaussian integration. 

In the non-zero-mode sector, relying on the oscillator algebra it is not difficult to 
compute^ 

(n.z.)(-Bop,0|e~ "* "|-Bop,-^)(n.z.) 



n=l 

oo 

--q 48 |j(l-g«) 2 =[i^{t)] 2 , 



n=l 

where we introduced, for later notational convenience, 

q = e-4-t = e2-i- (4.19) 

and wrote the result in terms of the Dedekind eta function (see Appendix A for our con- 
ventions). Recall that t is a real parameter. 

The matrix element in the ghost sector is such that it cancels exactly the non-zero- 
mode contribution of two directions. Thus, putting together all the pieces, we obtain 

W{L,R) = V2^\N\'' ^e-'^^'*-^ [r/(r)]-^ . (4.20) 

Jo f+2 

Let us consider the expression yV(i?, L) obtained exchanging L and R. Let us change 
integration variable in this expression by setting t = l/(4t'), which corresponds, according 
to eq. (4.19), to r = —1/t'. We can then use the modular properties of Dedekind's 
function, see eq. (A. 3), and find 

W(i?, L) = ^/2^|AA|^22--i-^ r ^^^ e-«^*'-^ W)] "^ . (4.21) 



^Notice that all the D + 2 direction contribute in the same way to the matrix element, despite the fact 
that the exponent in the boundary state 4.16 has the opposite sign in front of the {a]^)"^ oscillators with 
respect to the other ones. 
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This coincides with }V{L,R), as given in eq. (4.20), provided we fix 

D 1 , , 

. = - + -. (4.22) 

4.3 Explicit integration and loop expansion 

The integral over the parameter t in the definition of the Wilson loop amplitude can be 
explicitly performed, by techniques very similar to those exploited in [26, 27] to handle 
other geometries in the same operatorial approach to the effective string. We can Fourier 
expand the power of Dedekind's eta function appearing in eq. (4.21) writing 

OO OO CO 

ivir)]-^ = E Cke-'-<'-^s) =Y^c, ,^-il =Y,cuqK (4.23) 

fe=0 fc=0 k=0 

where for simplicity we introduced 

k = k-— . (4.24) 

48 ^ ^ 

Taking into account eq. (4.22) we have then 

W{L,R) = V2^\Ar\'f:c, r^t-fe--[^^^^*^^^] . (4.25) 

fc=o -^0 * 

We can now express the integral in terms of modified Bessel functions^ and write 

OO . s H 

W{L,R) = 2V2^\M\''^cJ-^] K^iaASk), (4.27) 

k=0 V ^ / 

where for brevity we set 

(y = 
2 

and introduced 



D 

a = - (4.28) 



£k = ]Jl + ^k. (4.29) 

Eq. (4.25) is well suited for an expansion for large area (measured in units of the 
inverse string tension), i.e., for large a A, keeping finite the ratio u = L/R. As we saw 
in section 2, in the physical gauge Nambu-Goto description this corresponds to the loop 
expansion of the sigma-model. Using the asymptotic expansion of the modified Bessel 
functions for large values of their arguments, see eq. (A. 13), and the expansion of £k we 



We use the formula 



t \B 



t-'c-'' '-— = 2 - K-,{2AB) . (4.26) 
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find 



WiA,u) = V2^\M\' 



2\ 4 



vr 



2aA 



e-'^'^ Y^ Cfce-2-«^ 



1 



X Sl + — r 

I aA 



-a^-4 a -2 



+ 



-TTuk — 2tt u k 



+ 



[aAy 



32 2 

a'' - 40q2 + 144 a^ - 60^ - 4a + 24 



2048 



+ 



64 



TTuk 



+ ^ ^^" + ^^ 7rVP - (a - 6)7r^u^P + 27r4«4fe4' 

16 



+ 



[aAY 



(4.30) 



Note that, recalling eq. (4.23), we have 



E 

k=0 



c,e-2-"^ = J]Q^ = [7?(in)]- 



(4.31) 



the introduction of Q = e"^'^" is convenient for later manipulations. The powers of k 
appearing in eq. (4.30) correspond then to powers of the logarithmic derivative Q{d/dQ) 
applied to the sum in eq. (4.31). In other words, we have 



W{A,u) = \/2^|AAp 



2\~ 



■K 



-a A 



2aA 



M'+^ 



a^ -A a -2 



32 



.d 



+ ^'"'l«5c3j+^''"'i«5Q 



,d 



+ 



[aAf 



a^-40a^ + U4 a^-Qa^-Aa + 2A (d 
mS + 64 "" V^dQ 



+ 



3(q2-8q + 12) 2 2 ( r^d 



16 



dQ 



^^ Q-J7^ + (a - 6)7r^n^ Q-jp: + 2ti^u^ Q 



.d 



dQ 



'-Y 

dQ) 



O 



{aAr 



[r](}u)\ 



(4.32) 
Logarithmic derivatives of the r/ function can be expressed in terms of Eisenstein series, see 
Appendix A for details. By doing so, recalling the expression eq. (4.28) of the parameter 
a in terms of the dimension D, we finally get 



W{A,u) = \/2^|AAp 



_^...(i)*,,.„„-f{.,^,|M,...},,.33, 



where 



C2{U) 



V24 



D{D - A) 



I — I I 2Du^Ei{iu) ''-^-^E2(iu)E2(l/u) ] + ^ ^^^- ^ (4.34) 



{D + A){D-A) 



128 



18 



and 



+ 



24/ 
D{D- 



{D + 2A)D 

6 

A){D-8){D 



u'El 



m 



3Z:>(D - 8)(Z? - l2)u^Ei{iu)E2{iu)E2{i/u) 



■Elim)Elii/u) 



+ (24 



D{D- VI) 



u-'Eg{iu) 1 



TT 



-^^2(1^^; 



—D(D 

L64 ^ 



(i:) + 4)i?(Z) 



256 



4)(L>- 12)u2E4(m) 



+ 



(D + 12){D + 4)(Z) - A){D - 12) 
32768 



(4.35) 
Using the transformation properties of the Dedekind function and of the Eisenstein 

series given in Appendix A it is not difficult to check that eq. (4.33) is invariant under the 

transformation u 1/u, i.e., under the exchange L -H^ R. 

The term in eq. (4.33) that multiphes the curly brackets should correspond to the 

classical plus one-loop result discussed in sec. (2). Up to numerical normalization, we have 



W(i_ioop) oc A 



'^u s [r]{m)y 



^-oA ^ ^-(i+f ) 



rj{\u) 



D_ 

2 



tA 



(4.36) 



By comparison^" with eq. (2.5), we deduce that the observable W(L,i?) of the first-order 
formulation is in fact proportional to the Wilson loop W{L^ R) through a modular-invariant 
prefactor factor of {LR)~'^ = [LR)^^2'^T>. ^From the operatorial treatment we extract 
thus the following prediction for the Wilson loop up to three loops: 



W{L,R) oce 



-a A 



ry(iu) 



R 



ID 

' 2 



1 + 



C2{u) C3{u) 



aA 



+ 



[aAY 



+ 



(4.37) 



The two-loop coefficient £2(1*) coincides with the (corrected) expression eq. (2.7) for the 
two- loop term C2, up to the constant term^^. This very remarkable and non-trivial agree- 
ment makes us confident in our result, and in particular in the third loop corrections 
contained in eq.s (4.37,4.35). Notice that the overall normalization of eq. (4.37) is totally 
irrelevant for the observable R^{L,R,n) which we have simulated numerically, since it is 
given by the ratio of two Wilson loops. 

We can also utilize the exact expression eq. (4.25) for W{L,R), together with the 
rescaling of the one-loop term just discussed, to argue that the Wilson loop can be written 
as 



W(L, R)ocA'^y"cki — ] Kd (a as, 



k=0 



£k 



(4.38) 



This expression can be used to estimate (by truncating the series) the all-loop prediction 
of the NG effective string for the Wilson loop and for the observable R^{L, R, n). 



Note that in this section we are considering the first order formulation of the NG action only, without 
boundary terms, and thus we do not obtain the perimeter term. 

^^In the comparison with the data, the eflfect of this constant is negligible in all cases except the one 
reported in fig. 2(d), where it enhances the agreement. It would be interesting to understand better this 
difference with respect to the Dietz and Filk treatment. 
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Figure 4: Plot of R^{L, —,n) for various values of n and n against the quantum string 
corrections up to three loop and the whole Nambu Goto prediction. 



5. Beyond two loops: comparison with the data 

Using eq.s (4.37,4.35) we may easily extract the three loop correction for the observable 
R^ in which we are interested. Similarly, using eq.(4.38) we may extract the correction 
which one would obtain assuming the validity of the Nambu-Goto action to all orders. 
These predictions are reported, together with the numerical data in fig. (4). To allow a 
simpler comparison we reported in fig. (5) an enlarged version of the plots, restricted only 
to the data for i? > 20 which, as discussed in section 3.3, are not affected by boundary 
corrections. 

Remarkably enough, the three loop correction turns out to have the opposite sign with 
respect to the two loop one and thus it goes exactly in the same direction as the numerical 
data. The statistical errors are too large to allow any reliable test but, as it can be seen 
in fig. (4) the gap between two and three loop corrections increases as R decreases and in 
principle, once the boundary corrections will be under control, this particular combination 
of Wilson loops could become a perfect setting for high precision tests of the universality 
conjectures. In this framework it is also interesting to notice that this change of sign is 
somehow peculiar of the three loop correction and is not present in the whole Nambu-Goto 
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Figure 5: Same as fig. (4, but keeping only the data for i? > 20 whicli are not affected 
by boundary conditions (see tlie main text). To make easier tlie comparison between the 
predictions at different orders, we report here the ratio of E1„ to the one-loop prediction; 
the error bars are rescaled accordingly. 



correction (which is also reported in fig. (4) for comparison). While a good agreement 
of the data with the three loop correction is certainly expected due to universality, the 
whole Nambu-Goto action should instead be excluded: it cannot be the correct effective 
action, by the arguments reviewed in the Introduction. The data seem indeed to support 
this statement, even if, as for the three loop case, the statistical errors are too large to 
allow any reliable statement and a definitive answer will require a precise control of the 
boundary correction in order to study the corrections in the range R < 20. 

6. Conclusions 

In this paper we have considered the disk partition function for the Nambu-Goto effective 
string theory corresponding to a rectangular Wilson loop. Our aim was to check the 
prediction of the NG model including its quantum corrections up to three loops against 
a numerical simulation tailored to this goal. This purpose, which goes beyond the tests 
available in the literature, owes its interest to the theorem which states that the corrections 
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up to three loops are universal for all effective string models; a check of the NG model to 
this order would in fact represent a stringent test of the very idea of an effective string 
description. 

Setting up such a test required progresses both on the theoretical side and on the side 
of simulations. 

On the theoretical side, we have extended to the Wilson loop geometry an operatorial 
approach based on the first order re-formulation of the NG model that had been already 
used for the cylinder and torus geometry. The operatorial evaluation makes use of open 
string analogues of boundary states and leads to an exact expression. This result resums 
the loop expansion that arises in the physical gauge approach, and can be expanded to the 
desired loop order. In particular, we have determined the three-loop correction, which had 
not been obtained in the physical gauge approach. 

On the numerical side, we set up the Montecarlo simulation of an observable corre- 
sponding to the ratio of two Wilson loops. This observable is devised so as to minimize 
the effects of boundary terms in the effective action, for which theoretical predictions are 
not available. Our simulation has been carried out in the 3d Ising gauge model, and has 
reached a level of precision such that it confirms nicely the validity of the effective string 
approach up to two loops; this already represents a big improvement with respect to the 
available literature. The statistical errors are yet too large to test reliably the three loop 
correction, even though it seems to go in the right direction. 
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A. Notation, conventions and useful formulae 

Here we establish our notations and give some properties of modular functions and Bessel 
functions that we use in the main text, expecially in section 4. 

Dedekind function and Eisenstein series Dedekind's ry-function is defined, in terms 
of the quantity q = exp{27rir}, by 



v{t) = g^ n (1 - 9") • (A-i) 



n=l 
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One can expand its inverse in g-series: 

oo 

[v{rT'=Y.Pkq''-^^ , (A.2) 

fc=0 

where pk denotes the number of partitions of k. Other powers of r] admit a similar Fourier 
expansion, and we make use of this fact in eq. (4.23). 

Under the generators of modular transformations the Dedekind eta function transforms 
in the following way: 

Vir + l) = ei^r?(r) , 

ry(-l/r) = (e~'2r)2r/(r) . 

Eisenstein series can be defined through their Fourier expansion, which takes the form 

oo 
^ ' n=l 

where o"2fc-i(ra) is the sum of the {2k — l)-th power of the divisors of n: 

(7«(n) = ^(i". (A.5) 

d\n 

In particular, we have thus 

oo 

£;2(r) = l-24^ai(n)g2«, 

n=l 

oo 

E4(r) = 1 + 240^ (73 (n)g2«, (A.6) 

n=l 
oo 

E6(r) = l-504^c75(n)g2" . 



ra=l 

For /c > 1, the Eisenstein series are modular forms of weights 2k: 



^2fc \c7Td) = ^^^ + DY''E2k{T) , (A.7) 

so that in particular 

i?4(-l/T) = T^E^ir) , i?6(-l/r) = T6i?4(r) . (A.8) 

The series E2{t), instead, is almost a modular form of degree 2, since its transformation 
under the 5-generator has a non-homogeneus term: 

E2{--) = r^E2ir)--iT. (A.9) 

r vr 

Taking into account these transformation properties, it is not difficult to check that the 
various expressions given in section 4.2, and in particular eq. (4.33), are invriant under the 
exchange L ■h^ R, i.e., u -H- 1/u. 



23 



The E2 series is related to Dedekind r/-function by 



—drlogr]{T) = qdqlogr]{T) = — 



—drlogr]{T) = qdglogviT) = —,E2{t) . (A.IO) 



Further derivatives connect £'2 to £'4 and -Eg- 



qdgE2 = -^{Ei-El) , 

qd.Ei = -^{Ee - E2Ei) , (A. 11) 

qdgEe = -^{Ej - E2Ee) . 

Applying these formulae, we can evaluate the multiple logarithmic derivatives with 
respect to Q = exp(27rin) that appear in eq. (4.32): 

Q^) [^M]"" = -^ [^M]""^2(m) , 
2 






•S^T^ ) [^M]"" = -77T7\3 [^M]"" ( 16£6(m) + 6(a - 4)£4M£2M 



+ (a2-6a + 8)E|(in 

^h) ^''^'''^^~" ^ W ^''^'''^^'" f^^^^" " ^)^6M£2M + 12(a + 12)Elim) 

+ 12(a^ - 10a + 24)^4 (in) E|(in) + (q^ - 12q2 + 44q - 48)E|(in) 

(A.12) 

Asymptotic expansion of Bessel Functions The asympotic expansion of the modified 
Bessel functions of the second kind Ka{z) for large values of their argument z takes the 

form 

vr .^■^ 1 T{a + m + I) 1 



Kjz)^J-e-^y^ ;'^'^^"'^']> \ . (A.13) 

^^mlTa-m + i 2z™ ^ ' 
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